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1. Introduction 

The years since the COBE observations of the CMB have ushered in an age of precision cos¬ 
mology. Key cosmological parameters have been determined by measurements of the distribution 
of matter in the Universe through WMAP and Planck observations of CMB anisotropies and large 
volume galaxy surveys such as SDSS. These surveys have made precision measurements of param¬ 
eters describing the matter content of the Universe - the baryons Q.h, dark matter Clc, dark energy 

radiation and neutrinos fly - and the physics of inflation - via the tilt amplitude As, 
running dn^/dlogk or the primordial potential power spectrum and r the ratio of tensor-to-scalar 
modes produced by inflation. These measurements have firmly established a working model of 
our Universe, known widely as the ACDM model of cosmology. The success of this model is de¬ 
spite our ignorance of the physics of key components, such as dark matter and dark energy, and 
deviations from the standard model could help refine our undersfanding. 

Despife fhis precision, measuring model paramefers is only fhe flrsf sfep towards a deep un- 
dersfanding of fhe underlying physics. Our ignorance of fhe nafure of fhe dark mailer and fhe dark 
energy or how neulrinos acquire mass and whal value lhal mass lakes are jusl Iwo queslions lhal 
modem cosmology hopes to address. Over fhe nexl decade Iwo palhs will help shed lighl on Ihis. 
The simplesl is simply to measure Ihese cosmological parameters ever more precisely and over 
a wider range of limes and scales in Ihe hope of gaining furlher insighls. The exemplar of Ihis 
is wilh dark energy, where allempls to measure Ihe redshifl evolulion of Ihe dark energy density, 
paramelerised by an equation of slate w{z), mighl distinguish a Irue cosmological conslanl from 
more general dark energy or modified gravity. For olhers Ihere are critical Ihresholds of precision 
required to distinguish physical scenarios - for example, measuring Ihe sum of Ihe neulrino masses 
My <0.1 would determine Ihe neulrino mass hierarchy. Clearly more precision is a good Ihing, bul 
if is nol Ihe only palh forward. 

More generally, we can seek signalures of new physics in ways dislincl from Ihe dislribulion 
of large scale mailer. For example, Ihe processes lhal produce dark mailer will also allow il to anni¬ 
hilate and maybe to decay. The associated release of energy could have impacl on Ihe surrounding 
environmenl, heating Ihe inlergalaclic medium. Pursuing unique signalures of new physics in new 
regimes will be a key part of Ihe nexl decade. The SKA is uniquely placed to probe cosmology in 
Ihis way, as il is capable of mapping Ihe Universe over wide volumes and an unprecedented range 
of redshifls (see Figure 1). In Ihis chapter, we will focus on Ihe new opporlunilies created by SKA 
observations of the epoch of reionization (EoR) and the cosmic dawn (CD). This period has never 
before been observed offering a unique opportunity to test the consistency of the ACDM model 
and search for new hints to the great unanswered questions of cosmology. 

Fundamental physics in cosmology is generally associated with the density field, whose fluc- 
lualions are generated by inflation and which conlains imprinl of olher physics such as neulrino 
mass. Aslrophysics is a major challenge to gelling al cosmology wilh SKA, bul we can identify 
several key approaches to exlracling cosmology: (1) direclly from density fluclualions (2) via Ihe 
presence of exotic sources of radiation (3) via Ihe radiation fields produced by all sources since 
Ihose sources will Irace Ihe density field in some biased fashion (4) via Ihe weak lensing of Ihe 
21cm signal by slruclures belween Ihe observer and signal (5) miscellaneous olher probes. Ex- 
Iracling cosmology from Ihe 21cm signal during Ihe EoR will require innovative new techniques 
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Figure 1: Illustration of the volume probed by SKA where 3D comoving volume has been mapped to a 2D 
disk. The volume probed by an all sky galaxy survey out to z = 0.3 (red circle, ~SDSS) and z = 2 (blue 
circle, ^SKA-MID) are marked as is the volume at the redshifts probed by SKA-LOW, z = 6 — 27 (cyan 
annulus). Existing cosmological parameters are derived using the CMB and relatively local galaxies, with 
an implicit assumption that nothing strange happens in between. 

to separate astrophysics from cosmology. Nonetheless the sensitivity of the instruments, large vol¬ 
ume probed, and new redshift regime accessible to SKA makes this a very interesting area for new 
science. 

In this chapter, we will explore these different avenues for extracting cosmology from the 
21cm signal and attempt to assess the sort of constraints that will be achievable by SKA Phase 1 
and 2. However, we caution the reader that this is not a settled area and it is still unclear how well 
astrophysics can be dealt with. New ideas may improve the constraints, but new obstacles may 
render them optimistic. 

2. Cosmological parameters from density fluctuations 

In this section, we explore the ability of SKA to constrain cosmological parameters via ob¬ 
servations of the density field. Just as galaxy surveys constrain cosmology by using galaxies as 
a tracer of the linear density field, SKA can consfrain cosmology by using fhe 21 cm brighfness 
femperafure as a fracer of fhe densify field. This is nol an unproblemafic asserfion, since brighfness 
femperafure flucfuafions may be sourced by variation in fhe spin femperafure and neufral fraction 
in addition fo fhe densify field. 



( 2 . 1 ) 


Equafion 2.1 shows how fhese differenl terms come info play (Furlaneffo ef al. 2006). In a regime 
where Tg S> Tcmb and xhi = 1 then 5Tb will be an unbiased tracer of the density field. Af all 
ofher fimes fhe effecfs of asfrophysics musf be modelled and removed or somehow avoided. One 
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possibility might be to exploit redshift space distortions that produce an angular dependence of the 
power spectrum, which in the simplest linear theory models look like 


P{k)=P^0{k)+P^2^^+P^,^\ 


( 2 . 2 ) 


In principle, measurement of this angular dependence of the power spectrum could separate cos¬ 
mology and astrophysics since the P^4 = Pg so directly probes the density field. In practice, this 
separation is complicated by non-linear growth of structure (Shaw & Lewis 2008; Mao et al. 2012) 
and the motion of ionised regions themselves (McQuinn et al. 2006) and it is unclear how effective 
it can be. We will return to a discussion of separating astrophysics and cosmology in §3 as this is a 
critical point. 

In this section, we take the optimistic view that there will a regime in which 5Th oc (1-1-5) so 
that the 21cm signal provides a clean measurement of the density field. This approach enables us 
fo evaluafe fhe besf case scenario for SKA in measuring cosmological parameters. By comparing 
fhis fo galaxy surveys we gel a sense of how compelilive SKA could be, if aslrophysics could 
be overcome. While we focus on sfandard cosmological parameters - energy densily in baryons 
D.h, mailer cosmological conslanl bubble parameter /i; inflalionary parameters Ag, jig and 
dus/dlogt, neutrino mass My and curvature - SKA will open a new regime into exotic physics 
that can only be probed at high redshift, for example compensated isocurvature modes whose effect 
decreases with time (Gordon & Pritchard 2009). Cosmology is moving from simply wanting to 
measure cosmological parameters more accurately and instead becoming more focused on control 
of systematics and relaxing simplifying assumptions. SKA will test consistency of cosmological 
parameters in a new redshift range. 

The sensitivity of a radio interferometer to the 21cm power spectrum has been well studied 
(e.g. Bowman et al. 2006; McQuinn et al. 2006; Mao et al. 2008; Mellema et al. 2013) and we 
follow the same approach here. The variance of a 21 cm power spectrum estimate for a single 
k-mode with line of sight component kn = /ik is given by (Lidz et al. 2008): 



(2.3) 


The first term on the right-hand-side of the above expression provides the contribution from 
sample variance, while the second describes the thermal noise of the radio telescope. The thermal 
noise depends upon the system temperature Tsys, the survey bandwidth B, the total observing time 
tint, the conformal distance D{z) to the center of the survey at redshift z, the depth of the survey AD, 
the observed wavelength A, and the effective collecting area of each antennae tile Ag. The effect of 
the configuration of the antennae is encoded in the number density of baselines n±{k) that observe 
a mode with transverse wavenumber k± (McQuinn et al. 2006). Observing a number of fields A^fieid 
further reduces the variance. Given the sensitivity of the instrument, the Fisher matrix formalism 
can be used to estimate 1 — a errors on the model parameter A; are (Fr.^)'/^, where 



(2.4) 
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In this equation, ^survey = D'^AD{X'^/Ag) denotes the effective survey volume of our radio tele¬ 
scopes and we assume wavenumber bins of width Ak = ek. Key to determining cosmological pa¬ 
rameters are the effective volume probed and the minimum wavenumber probed ^min where modes 
can still be assumed to be linear. SKA has a significant advantage over galaxy surveys as more 
modes are still in the linear regime at z > 6. 

Table 1: Low-frequency radio telescopes and their parameters. We specify the number of antennae Na, total 
collecting area Atot, bandwidth B, and total integration time tint for each instrument. These values are fixed 
at V = 1 lOMHz and extrapolated to other frequencies using Atot = iVamiVdip^dip with a physical station size 
of 35m and the number of antennae per station Nup = 289 and A^ip = min(A^/3,3.2m^). 


Array 

Na 

Atot (10^ m^) 

B (MHz) 

tint (hr) 


^max {klfl) 

MWA 

112 

1.6 

8 

1000 

4 

0.75 

PAPER 

128 

0.9 

8 

1000 

4 

0.15 

EOEAR Core 

48 

38.6 

8 

1000 

100 

1.5 

HERA 

331 

50.0 

8 

1000 

14.3 

0.3 

SKAO 

899x0.5 

831x0.5 

8 

1000 

35 

2 

SKAl 

899 

831 

8 

1000 

35 

2 

SKA2 

899x4 

831x4 

8 

1000 

35 

2 


We first illustrate the sensitivity of different iterations of SKA in Figure 2, where we take the 
parameters in Table 1 for SKAO - with 50% of the SKAl baseline collecting area, SKAl, and 
SKA2 - with x4 the collecting area of SKAl. For each of these we assume a filled core followed 
by distribution out to a maximum radius /?max- HERA is assumed to have a uniform antennae 
distribution. SKAl has 911 stations total with 899 in the core and 650 stations within a radius of 
1km accounting for ~75% of the total number of stations and collecting area (Dewdney 2013). At 
lower frequencies the array is densely packed and has constant collecting area, at higher frequencies 
the array becomes sparse. 

Figure 2 illustrates a few key points governing parameter constraints. Here we have elim¬ 
inated modes whose wavelength exceeds the instrument bandwidth removing sensitivity to the 
largest physical scales (smallest k modes). At z = 8, SKAO is directly comparable in sensitivity 
to the proposed HERA experiment (Pober et al. 2014), which is more centrally concentrated to 
compensate for its small number of stations. Detection of the 21cm signal at z > 20 with SKAl 
is dependent upon either a strong 21cm absorption signal that boosts the amplitude of the 21cm 
power spectrum, e.g. Tg <C Ty as expected before X-ray heating, or spin temperature fluctuations 
that add additional power over that of the density field. Unfortunately, it seems likely that dur¬ 
ing the absorption regime the details and spatial variation of the spin temperature will matter and 
complicated getting at cosmology. 

Table 2 shows the cosmological parameters obtained with the listed experimental perfor¬ 
mances using a Eisher matrix approach following McQuinn et al. (2006). The key take home 
message of this is that SKA-EOW has the raw sensitivity to add useful information on cosmolog¬ 
ical parameters to Planck. The largest gains are on parameters that require small scale informa¬ 
tion, for example the running of the primordial power spectrum and related inflationary parameters 
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Figure 2: Sensitivity plots of HERA (red dashed curve), SKAO (red), SKAl (blue), and SKA2 (green). 
Dotted curve shows the predicted 21cm signal from the density field alone assuming jt// = 1 and Tg Tcmb- 
At z = 20, we also plot the case of Tg = 20K in the z = 20 panel to give a better sense of the expected 21 
cm signal during absorption. Vertical black dashed line indicates the smallest wavenumber probed in the 
frequency direction k = 2n/y, which may limit foreground removal. Left panel: z = 8 Right panel: z = 20. 


(Barger et al. 2009; Adshead et al. 2011) and neutrino masses (Pritchard & Pierpaoli 2008). This 
also indicates that SKA-LOW will have the sensitivity to provide a useful consistency check on 
cosmological parameters from the high redshift regime long before dark energy becomes impor¬ 
tant. 

These numbers assume a single deep field designed to reduce thermal noise and so maximise 
sensitivity on the smallest scales. This tends to maximise the constraint on parameters like neutrino 
mass, which modify the power spectrum primarily on small scales. On large scales, cosmic vari¬ 
ance dominates over thermal noise. This makes it useful to complement a single deep field with 
many shallower fields, which increase the survey volume and reduce the cosmic variance. The 
SKA-LOW survey strategy of shallow ~ 10000 deg^, mid 1000 deg^, and deep 100 deg^ surveys 
provides a good mix to optimise for cosmology. 

Table 2: Fiducial parameter values and 1 — a constraints on cosmological parameters. Non-cosmological 
parameters included in the analysis {T, x//(z = 7), v//(z = 7.5), v/r(z = 8)} are not shown. We take = 
2Mpc^* as the limit to linear modes. 


Value 

\ogCl,„h^ 

-1.9 

logClhh^ 

-3.8 

CIa 

0.7 

tls 

0.95 

log(A^/10 1°) 
-0.19 

Clk 

0 

dris/dlogk 

0 

Mv (eV) 
0.3 

Planck 

0.028 

0.0068 

0.038 

0.0035 

0.0097 

0.0022 

0.0047 

0.35 

Hera 

0.0091 

0.0055 

0.011 

0.003 

0.0088 

0.0021 

0.0036 

0.12 

SKAO 

0.017 

0.0058 

0.023 

0.0032 

0.009 

0.0022 

0.0034 

0.22 

SKAl 

0.0083 

0.0051 

0.01 

0.003 

0.0084 

0.002 

0.0018 

0.12 

SKA2 

0.0016 

0.0048 

0.0026 

0.0027 

0.0081 

0.0012 

0.00092 

0.084 


We make no attempt here to model the effects of astrophysics on these constraints. Increas- 
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ingly conservative assumptions can degrade these constraints arbitrarily far (Mao et al. 2008), so 
these should be viewed as optimistic bounds on the constraints that might be achieved. Nonethe¬ 
less it is clear that the attempt to extract cosmology from CD and EoR observations could be quite 
rewarding. 

3. Separating “gastrophysics" and cosmology 

The key challenge for extracting fundamental physics from the 21cm signal will be separat¬ 
ing the effects of cosmology from “gastrophysics". A number of avenues have been studied in 
the literature, which broadly separate into (1) avoidance and (2) modelling and (3) redshift space 
distortions. The optimistic case in the previous section assumed the possibility of avoidance - a 
clean region in redshift where Tg 3> Tcmb and x// = 1. Theoretical modelling of the 21cm signal 
(e.g. Pritchard & Loeb 2008; Thomas & Zaroubi 2011; Mesinger et al. 2013; Fialkov et al. 2014) 
suggests that we are unlikely to find such a region, although certain epochs may approach this 
limit. In the absence of a clean window, it might still be possible to avoid astrophysics via the 
angular dependence of the power spectrum induced by redshift space distortions. Focussing on the 
P^4 Ri Ps part could lead to clean cosmological measurements. Obtaining precision cosmology this 
way is hard and the literature suggests little improvement over Planck will be possible (McQuinn 
et al. 2006; Mao et al. 2008). Figure 3 shows predicted errors bars for SKA on the P ^2 and P^a 
parts of the power spectrum. A detection is possible with at wavenumbers k = 0.1 — 1 Mpc^^, but 
with much less precision than the full 21 cm power spectrum. 



Figure 3: Sensitivity plots at z = 8 on P ^2 (top panel) and P^a (bottom panel) for HERA (red dashed 
curve), SKAO (red), SKAl (blue), and SKA2 (green). Dotted curve shows the predicted 21cm signal from 
the density field alone assuming xh = I and Ts ^ Tcmb- Vertical black dashed line indicates the smallest 
wavenumber probed in the frequency direction k = 2n/y, which may limit foreground removal. 


The most likely path is to model the contribution of astrophysics. Compared with the CMB 
our theoretical understanding of the 21cm signal during reionization is poor. Predictions for the 
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21cm power spectrum do not exist at the same level of precision as the cosmology. Nonetheless, we 
expect the contribution of astrophysics to be relatively broad band and determined by extra power 
about a characteristic scale, eg the bubble size during reionization. Mao et al. (2008) showed that 
relatively simple parametrisations capture the shape of ionisation contributions and so might be 
fitted for and marginalised out. Information from measurements of would complement this, 
as would information from different redshift slices. Given the large amount of information in 3D 
and the ability of SKA to image the signal - allowing ionised bubbles to be directly identified and 
masked ouf - if may be possible fo characferise fhe asfrophysics on large scales. This has yef fo be 
examined in defail and if is unclear how far 21cm observafions mighf be “cleaned" of asfrophysics. 

One fhing fo nofe is fhaf reionizafion desfroys information - ionised bubbles produce no 21cm 
signal - while healing and Lya coupling merely overlay fhe densily field wilh olher information. If 
may be possible in fhe fulure fo eslablish a way of separaling spin lemperalure fiuclualions from fhe 
densily field in some olher way, e.g. by using galaxy observafions fo reconslrucf fhe Lya flux on 
large scales, and so recover more of fhe cosmological informalion from observafions in fhe cosmic 
dawn. 

4. Constraining new physics from heating 

The 21cm signal probes bolh fhe ionizalion and Ihermal slafe of fhe IGM. Allhough we do nol 
know fhe precise liming and evolution of fhe signal, empirical scaling relalions based on local slar- 
forming galaxies (e.g. Mineo el al. 2012) suggesl lhal fhe X-rays from early galaxies heal fhe IGM 
fo femperafures above fhe CMB before fhe bulk of reionizafion (e.g. Furlanello 2006; McQuinn & 
O’Leary 2012). This marks fhe Iransilion of fhe 21cm signal from absorption fo emission, wilh 
large-scale flucfualions in gas lemperalure likely driving fhe 21cm power fo ifs largesl amplilude 
(e.g. Prilchard & Furlanello 2007; Back el al. 2010). The epoch of IGM heating is a powerful 
probe of Ihe high-energy processes in Ihe early Universe, wilh could have bolh aslrophysical and 
cosmological origins. Bolh can fell us aboul Ihe nalure of dark mailer (DM). 

In order lo explain Ihe apparenl deficiencies of CDM on small (sub-Mpc) scales. Warm Dark 
Mailer (WDM) models have recenlly gained in popularity. In Ihese models, DM is assumed lo 
consisl of smaller mass particles, ~ keV, such as Ihe sterile neulrino or gravilino. The increased 
particle free-slreaming and velocity dispersion (acting as a sorl of effective pressure), can dra¬ 
matically suppress slruclures on small-scales. This suppression is even more obvious in Ihe early 
Universe, where typical halos hosting galaxies were much smaller, and larger slruclures did nol 
have time lo fragment Currenl aslrophysical lower limils on Ihe WDM particle range from nix > 
1-3 keV (assuming a Ihermal relic relativistic al decoupling), wilh various degrees of aslrophysical 
degeneracy (e.g. de Souza el al. 2013; Kang el al. 2013; Pacucci el al. 2014; Viel el al. 2013) 

The resulting dearlh of galaxies in Ihe early Universe means lhal Ihe aslrophysical epochs in Ihe 
21cm signal were delayed. The challenge as always will be lo disenlangle Ihe cosmological impacl 
from aslrophysical uncerlainlies, for example a lower lhan expected slar formation efficiency in 
CDM would look superficially similar lo a higher slar formation efficiency in WDM. Since Ihe 
fractional suppression of slruclure increases wilh redshift. Ibis becomes much easier wilh Ihe firsl 
galaxies observable wilh Ihe SKA. For example, we only need lo undersland Ihe asfrophysics of 
Ihe firsl galaxies lo an order of magnilude in order lo improve on currenl nix conslrainls (Silwell 
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et al. 2014). Moreover, even if the star-formation efficiency in CDM is allowed to vary in order to 
mimic the mean 21cm evolution in WDM models, the signal will still not be completely degenerate 
(see Fig. 4a). This is due to the fact that the galaxies driving the 21cm evolution in WDM should 
reside in higher mass, more rapidly evolving halos, than those in CDM. The increased bias of such 
halos results in a larger 21cm fluctuations (see Fig. 4a). 

The heating of the IGM could also have a cosmological component. In particular, annihilations 
of dark matter particles in the ~ 10 GeV mass range (motivated by recent results from indirect 
experiments; (e.g. Adriani et al. 2009; Abdo et al. 2010; Aguilar et al. 2013) could provide a 
dominant source of heat, before the birth of the first galaxies. Driven by the evolution of Mq 
structures, several orders of magnitude smaller than those hosting galaxies, heating is expected to 
be much slower in such models, resulting in a smaller brightness temperature gradient dST^/dv ~ 
4mKMHz^^ in the range V ~ 60 — 80 MHz (Valdes et al. 2013). Moreover, DM annihilations 
would heat the IGM quite uniformly, which is not the case for heating driven by astrophysical 
sources residing in early galaxies. The resulting lack of temperature fluctuations (see Fig. 4b) 
would result in dramatic drop in 21cm power during heating, which would be easy to identify with 
the SKA (Evoli et al. 2014). Furthermore, the ensuing rise in 21cm power when the galaxies start 
contributing to heating the IGM should occur when the IGM is already in emission. The later is a 
qualitatively robust signature of DM annihilation heating, easily obtainable with the SKA. 


CDM 

1.0 
0.8 
S 0.6 

u 

0.4 

0.2 
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Figure 4: Right panel: Evolution of the power spectmm of 57), for WDM with mx = 2 keV. The top panels 
show power spectra &tk — 0.08, 0.18 Mpc^* for WDM (dashed) and the CDM model (solid). CDM models 
have /*(z) (star-formation efficiency) chosen to reproduce the global 21-cm signal found for the respective 
WDM model. The bottom panels show the difference in the power spectrum between WDM and CDM 
models. Dotted curves show forecasts for the 21cm power spectrum thermal noise as computed in Mesinger 
et al. (2014) assuming 2000 h of observation time. The dotted green, blue and red curves are the forecasts 
for the MWA, SKA and HERA, respectively. This figure is from Sitwell et al. (2014). CDEs of Ty/Ts 
corresponding to the fiducial and extreme astrophysical X-ray heating (black and gray curves respectively) 
from Mesinger et al. (2013). The colored curves correspond to models in which 10 GeV DM annihilations 
are also accounted for (in addition to fiducial astrophysical heating), with varying relative contribution. The 
curves correspond to the redshift for which Ts ~ Tcmb- Eigure is from (Evoli et al. 2014). 
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4.1 21 cm signal from primordial magnetic fields 

Primordial magnetic fields (PMFs) has been intensively investigated in the literature as possi¬ 
ble seeds for large scale magnetic fields observed in galaxies and clusters of galaxies (for a recent 
review, see (Durrer & Neronov 2013)). Magnetic fields in galaxies in high redshifts (Bemet et al. 
2008) and in void regions (Neronov & Vovk 2010; Ando & Kusenko 2010; Takahashi et al. 2013) 
can well be the pieces of evidence that the seed fields are of primordial origin. The primordial mag¬ 
netic fields may be created in the very early universe, e.g., at the epoch of inflation, cosmological 
phase transition, and cosmological recombination. The Planck collaboration recently placed limits 
on the magnetic field strength smoothed on 1 Mpc scales Ba=impc < 3.4 nG and the spectral index 
ub <0 of any PMFs from the temperature anisotropies on large and small angular scales (Planck 
Collaboration 2013). 

The CMB brightness temperature fluctuations produced by the neutral hydrogen 21-cm line 
(21 cm) would offer a new probe of the primordial magnetic fields (PMFs) created in the early 
universe. For the 21 cm observation, aside from the early structure formation effect by the Lorentz 
force from the PMFs, one of the important effects is the dissipation process of the PMFs that in¬ 
creases the baryon temperature. The dissipation occurs mainly through ambipolar diffusion due to 
the velocity difference between neutral hydrogen (which is the dominant component in the dark 
ages) and ionized particles (whose trajectory is bent by the Lorentz force). The effect of the dissi¬ 
pation is rather significant. The gas temperature can reach 1000 K or even 10^ K at z = 30 if the 
magnetic fields have the strength of ~ 3 nG (Sethi & Subramanian 2005; Tashiro & Sugiyama 
2006; Schleicher et al. 2009; Kunze & Komatsu 2014). 

This dissipation will give rise to a unique signature of the PMFs on the 21 cm observation. 
Because the spin temperature is closely coupled to the gas temperature at high redshift (z > 30), 
the 21 cm signal would come as ‘emission’ if the energy dissipation is efficient. In Fig. 5 the global 
HI signal with several magnetic field strengths are shown. For cases with sufficient magnetic fields, 
say B > 0.03 nG, the signal is always emission against CMB while in the standard ACDM model 
the signal would be absorption for the frequency range of /v < 80 MHz (corresponding to the signal 
from redshift z > 20). 

We show the angular power spectrum of the 21 cm brightness temperature including the PMFs 
in Fig. 5b (Shiraishi et al. 2014). Here we do not account for any (standard) heating effects (i.e., 
UV, X, and La background emissions) to isolate and clarify the effects from the PMFs. On large 
scales which may be relevant to SKA observations, there are two distinct contributions. One is 
from the standard (adiabatic) density fluctuations enhanced by the heating from the PMFs, and 
the other is from the PMF induced density fluctuations dominant on smaller scales (Tashiro & 
Sugiyama 2006; Schleicher et al. 2009). We can see from the figure that B = I nG magnetic fields 
are marginally within reach for a statistical detection of the power spectrum. Stacking observing 
channels in principle will add more statistical power. 

The angular correlation function in real space including the effects from the PMFs is also 
studied in Sethi & Subramanian (2009). The function exhibits a distinct feature because the PMFs 
induce early structure formation and the small scale halos form more compared to the case in 
the standard ACDM model The signal from primordial magnetic fields shows oscillatory feature 
contrary to that in the standard ACDM since the matter power spectrum induced by the PMFs is 
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V [MHz] i 

Figure 5: Left panel: The global 21 cm signal with magnetic field strength B = 1 0.1, 0.03, and 0.01 nG 
(colored lines from top to bottom). The solid line corresponds to the model without primordial magnetic 
fields. Note that any other heating source than magnetic fields is neglected in the figure. Right panel: 
Angular power spectra for PMF strengths: B = 0, 0.1, 1.0 nG at z = 15. The bottom curve shows the 
power spectrum from the standard density perturbations for fully neutral medium without any heating and 
reionization processes. The red and blue curves correspond to the cases with heating by the PMFs with 
B = 0.1 nG and B = 1.0 nG, respectively. The heating induces deviations of the spin temperature from the 
CMB temperature and the signal is enhanced. The noise curves for SKAl and SKA2 with IMHz bandwidth 
are also shown as indicated. By courtesy of M. Shiraishi & H. Tashiro. For reference, at z = 15, A: « 
0.2(//1000)Mpc-'. 


blue and most of halos are formed at the scale close to the magnetic Jeans’ length. It has been 
argued that 5 sigma detection of the 0.5 nG magnetic fields will be possible with less than one 
week integration of SKA observation (Sethi & Subramanian 2009). 

5. Bulk flows 

Tseliakhovich & Hirata (2010) demonstrated the existence of coherent supersonic velocity 
flows between baryons and dark matter after decoupling at z ss 1100. This has the consequence of 
inhibiting the formation of star forming galaxies in low-mass halos (M < 10^ Mq) (e.g. Maio et al. 
2011; Stacy et al. 2011). If there is significant star formation in such halos, dependent upon H 2 
cooling and the absence of Lyman-Wemer background, then the radiation from such galaxies can 
lead to a significantly enhanced 21cm signal (Visbal et al. 2012; Fialkov et al. 2014; McQuinn & 
O’Leary 2012). The effect of bulk flows could then have a major effect on the 21cm signal during 
the period of Lyman alpha coupling and IGM heating. 

Although the details of this are still quite uncertain, if this enhancement exists, it opens the 
possibility of measuring cosmology at z = 20 — 27 with SKA. The relative velocity flows trace the 
sound horizon and so especially enhance the baryon acoustic oscillation feature in the 21cm power 
spectrum (see e.g. McQuinn & O’Leary 2012). The BAO signature provides a standard ruler, 
calibratable with CMB observations, to form an inverse distance ladder stretching from z = 1100 
through z ~ 20. Such measurements would strongly constrain the parameter space for departures 
from ACDM, such as early dark energy models (e.g. Bartelmann et al. 2006; Doran & Robbers 
2006). Imaging the 21cm structures induced by these bulk flows will be possible with both SKA- 
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LOW Phase 1 and SKA2. This will be a truly novel probe of cosmology at high redshift vastly the 
reducing the possibility that non-standard dynamics could hide in the observational void between 
low redshift galaxy surveys and the CMB. 


6. Cosmology on ultra-large scales with SKA-Low 

The measurement of very large scales provides an unique way to probe modifications to the 
standard cosmological model. In particular it is on scales past the matter-radiation equality peak 
that General Relativistic corrections become important, at scales above ~ 5 comoving Gpc/h (Jeong 
et al. 2012). Probing this region would allow to check for any inconsistencies in General Relativity. 
However, current surveys are still far from probing this region - as an example, the BOSS survey 
only probes scales up to 200 Mpc/h (Anderson et al. 2012). In fact, the total volume contained at 
z < 1 is sufficient to probe only ~ 8 modes with k = 27r/(5Gpc//j) while at 6 < z < 30 almost 200 
modes could be constrained reducing cosmic variance. 

Moreover, primordial non-Gaussianity will affect the clustering of biased tracers of dark mat¬ 
ter, by adding an extra correction Abx{z,k) to the Gaussian large-scale bias of a given biased 
tracer A: Abx{z,k) = 3[b^{z) — l]D.mHQ5c/[c^k^T{k)D^{z)]fni- Here, D.,n = Clb + ^dm is the total 
(baryons plus dark matter) matter fraction, Hq is the Hubble constant, 5c ^ 1.686 is the critical 
collapse density contrast of matter, T (k) is the matter transfer function versus the physical wave 
number k, and D+ (z) is the linear growth factor of density perturbations. Attempts at detecting this 
effect with redshift surveys have led to some constraints on /„/ (Giannantonio et al. 2014). 

By observing at very high redshifts, a low frequency interferometer would be able to probe 
these large 3-dimensional scales during the Epoch of Reionization and beyond. Although astro¬ 
physics will generate model dependent features on the 21cm power spectrum, it is expected that 
on large enough scales the power spectrum should follow the dark matter one with a different am¬ 
plitude. For instance, during the epoch of reionization, the ionisation power spectrum should be a 
biased linear tracer of the dark matter one on scales much larger than the bubble size (see Figure 6. 
Although the bias itself will depend on the assumed astrophysical model, measurements on these 
scales will allow to pick any scale dependence generated by primordial non-Gaussianity or General 
Relativistic corrections. 

By accessing these large volumes, probes of the high-z 21cm signal will not only measure 
these large scales but also have enough modes to reduce cosmic variance on the scales of interest. 
Having large fields of view will therefore be a key factor for these experiments. The field of view 
of SKAl-Fow ranges from 7 deg^ al 220 MHz (z ~ 5.5) fo 133 deg^ af 50 MHz (z ~ 27), going 
basically as (1 +zY- These correspond fo scales of ~ 380 comoving Mpc al z ~ 5.5 fo ~ 2.3 
comoving Gpc al z ~ 27. As a further example, a fixed redshifl bin of 0.1 would evolve on Ihe 
same range from 50 Mpc fo 5.5 Mpc (decreasing wilh z). We see Iherefore lhal, allhough Ihe 
volumes accessible by SKAl-Fow are quile large, especially al very high-z, we would require a 
telescope wilh aboul 100 deg^ al z ~ 8 fo probe Gpc scales. Somelhing lhal will probably have lo 
wail for SKA2. Figure 7 shows Ihe conslrainls on Ihe primordial non-Gaussianily parameter for 
differenl telescopes. 
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Figure 6: Ionization power spectra with non-Gaussianity of the local form from numerical simulations. We 
show f„i = (0,20,100) (dot-dashed, dashed, solid) for efficiency = (5.8,3.0) (thin black, thick red) at 
z = 7.5, where xhi = (0.50,0.75). Analytical fits are in dotted lines. 



X jV>/ZlX 



Figure 7: Top'. Marginalized /„/ constraints for cases with noise (thick) and without noise (thin), which 
overlap for Omniscope. We consider a bandwidth of 6 MHz, but assume foregrounds can be removed 
on scales larger than k|| = InjiyB'). Bottom'. Marginalized /„/ constraints as function of bandwidth and 
number of antennae. The bandwidth limits the number of modes and largest scale probed along the LOS 
(via the survey volume V B and k™" 1/fi), whereas a larger number of antennae for fixed array density 

increases the survey resolution and number of perpendicular modes (via n{u±), on large scales Aant, and 
V-^ant)- The color coding is the same as for the top panel. 

7. Cosmic shear and the EoR 

It is possible that the EoR signal could be used to measure weak gravitational lensing. In 
Zahn & Zaldarriaga (2006) and Metcalf & White (2009) it was shown that if the EoR is at redshift 
z ~ 8 or later, a large radio telescope such as the SKA could measure the lensing convergence 
power spectrum. However a very large /sky and a very compact low frequency array was assumed 
by those authors. Here the calculation is repeated with parameters that are more consistent with 
current SKA baseline design. 

The current plans for a 25 square degree survey with SKAl-Eow will preclude making com¬ 
petitive measurements of the cosmological parameters through their effects on the weak lensing 
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power-spectrum because of sample variance (this is not true of the SKAl-Mid at lower redshift 
where the survey area will be much larger). It still might be possible to map the leasing convergence 
within the 25 square degree EoR survey area and a mix of wider-shallower and narrower-deeper 
observing modes can further optimise for weak lensing. This would allow us to actually “see” 
the distribution of dark matter in a typical region of the sky, something that is only possible with 
galaxy lensing around very atypical, large galaxy clusters. This would provide a great opportunity 
to correlate visible objects with mass and test the dark matter paradigm. 

The previously mentioned authors extended the Fourier-space quadratic estimator technique, 
which was first developed in Hu (2001) for CMB lensing observations to three dimensional ob¬ 
servables, i.e. the 21 cm intensity field The convergence estimator and the corresponding 

variance on the lensing reconstruction are calculated assuming that the temperature (brightness) 
distribution is Gaussian. This will not be strictly true during the EoR, but serves as a reasonable 
approximation for these purposes. Note that the lensing reconstruction noise contains the thermal 
noise of the telescope which -assuming a uniform telescope distribution- is calculated using the 
formula 


_ 


{271)X 


3)t2 
sys 


Btohsf} 


cover^max V 


(7.1) 


where the system temperature Tgys at high redshifts is dominated by galactic synchrotron radiation 
and can be approximated by Tgys = 60 x (v/300MHz)^^-^^K (Dewdney 2013), B is the chosen 
frequency window, Gbs the total observation time, Dtei the diameter (maximum baseline) of the core 
array, ^max(A) = 2%D[^\jX is the highest multipole that can be measured by the array at frequency 
V (wavelength A), and /cover is the total collecting area of the core array Acoii divided by 7r(Dtei/2)^. 

The advantage of 21cm lensing is that one is able to combine information from multiple red- 
shift slices. In Fourier space, the temperature fluctuations are divided into perpendicular to the line 
of sight wave vectors kj^ = 1/r, with r the angular diameter distance to the source redshift, and a 
discretized version of the parallel wave vector k\\ = ^ j, where .if is the depth of the observed vol¬ 
ume. Considering modes with different j independent, an optimal estimator can be found by com¬ 
bining the individual estimators for different j modes without mixing them. The three-dimensional 
lensing reconstruction noise is then found to be (Zahn & Zaldarriaga 2006) 


N{L,v) 


f dH [l-LQj + L-(L-l)C|,_i|/ 

h ^ (2;t)2 2C“Cf°\|,y 


(7.2) 


Here, C“ = Qj- +C^, where Cij = [T{z)YPij with T{z) the mean observed brightness tempera¬ 
ture at redshift z due to the average HI density and Pij the underlying dark matter power spectrum 
(Zahn & Zaldarriaga 2006). For SKAl-Fow we can consider a 1,000 hr observation time and we 
choose B = 8 MHz and ~ 40, but with multiple bands v that can be stacked to reduce the noise 
sothatAL = l/^[A(L,v)]-'. 

V 

At redshift ~ 8, we can assume the SKAl-Fow Baseline Design (Dewdney 2013) parameters 
of Acoii ~ 0.3km^ with maximum baseline Dtei = 4km, while for SKA we can consider Acoii — 
1.2km^. The estimated lensing noise is shown in Figure 8 along with the estimated signal. Here 
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is the convergence field power spectrum at = 8 and Nt the lensing reconstruction noise 
assuming a reionization fraction /hi = 1. 



Figure 8: The lensing convergence field power spectrum, C/'^, for sources at z = 8 is shown as a solid black 
line and lensing reconstruction noise A/ as dashed lines. The blue dashed line is for SKAl-Low with 10 
8 MHz frequency bins around z = 8 spanning the redshift range z — 6.5 — 11. The red dashed line is for 
SKA2 and the same frequency bins. The vertical line is approximately the lowest L accessible with a 5-by-5 
degree field. Where the noise curves are below C^’^, typical fluctuations in the lensing deflection should be 
recoverable in a map. 

These results show that it might be possible to map the lensing signal over a range of angular 
scales. This measurement greatly benefits from the larger collecting area that will come with SKA2 
(we also note that considering a more compact array, i.e. smaller Dtei, also improves the signal-to- 
noise). Note that the sensitivity of SKAl-LOW is such that larger area surveys at lower integration 
time do not significantly improve this picture. Although multi-beaming to get larger area at fixed 
integration time would help with cosmic variance on large scales. For SKA-LOW an optimal 
observing strategy includes a mix of shallow-wide surveys to beat down cosmic variance on large 
scales and deep-narrow surveys to measure small scales. The weak lensing power spectrum can 
be better measured for redshifts after reionization using SKA-Mid and the same 21 cm intensity 
mapping technique discussed, but over a much larger area of sky (Pourtsidou & Metcalf 2014). 

8. Conclusions 

SKA-LOW will provide the first window onto cosmological information from the epoch of 
reionization and cosmic dawn at z = 6 — 27. This makes it unique among the diverse array of 
future cosmological experiments, which are typically limited to redshifts z < 3. This new view 
of the Universe will test the standard cosmological model deep in the matter dominated regime, 
where, in principle, we believe we know the evolution of the Universe. The absence of deviations 
from ACDM would observationally confirm our current assumption that a single cosmological 
model applies from the CMB to the present. Deviations from ACDM predictions would signal new 
physics and provide a new way of learning about the Universe. 

Precision cosmology with SKA-LOW has potential to be very interesting, but is made difficult 
by astrophysical contamination. The large volume and number of linear modes of the density field 


15 






EoR/CD Cosmology 


Jonathan Pritchard 


accessible at high redshift will one day lead to a revolution in precision cosmological constraints. 
SKA-LOW Phase 1 will take the first step in this direction and a mix of shallow-wide and deep- 
narrow observing fields can help opfimise for cosmology; for example mafching a narrow-deep 
lOOOhr field over a single field wifh broad-shallow lOhr fields. The greafer sensifivify of SKA2 
will push significanlly furlher. The key limifafion will be our abilify fo separafe asfrophysics fhe 
density field fo measure cosmological paramefers. If is clear fhaf fhe efforf should be made, buf 
our existing undersfanding of fhe asfrophysics is still loo crude fo make robusl prediclions. SKA 
will confribufe fo more precise measuremenls of ACDM paramefers, buf crucially will fesf fhe 
parameters in a new regime. 

Much clearer is fhe abilify of bolh SKA-Low Phase 1 and 2 fo measure fhe Ihermal hislory 
of fhe Universe al redshifls z = 6 — 27. Never before measured, fhe Ihermal hislory conlains in- 
formalion aboul exotic physics - dark mailer physics, primordial magnelic fields, and more. The 
time dependence and spatial varialion of such healing is qualifalively differenl from lhal of galaxies 
providing a clear observalional signalure accessible fo SKA. Al fhe same lime, SKA could mea¬ 
sure fhe effecl of supersonic relalive velocily flows belween baryons and dark mailer. This probes 
much fhe same physics as fhe CMB and mighl allow BAO measuremenls of fhe sound horizon as a 
slandard ruler fo measure fhe geomefry of fhe Universe. 

As SKAl becomes SKA2, greafer sensifivify will allow wider sky surveys fo fixed deplh 
for a given infegralion time. These large volume observalions will access super-horizon physics 
and, wifh fhe SKA2, if will be possible fo search for relativistic correclions and primordial non- 
Gaussianily on > Gpc scales. Wider sky area and greafer angular resolution also improves fhe 
ability of SKA fo measure weak tensing of fhe 21cm background. Weak tensing will map fhe dark 
mailer field in represenlalive palches of fhe Universe on scales from degrees fo lens of arc minules 
providing a unique view of fhe growlh of dark mailer info fhe cosmic web. 

In ihis chapler, we have focussed on fhe main palhs fo cosmology wifh SKA-LOW. Many 
more speculalive ideas for cosmology exisl. For example, conslraining lime evolution of fhe fine 
slruclure conslanl (Khalri & Wandell 2007) or observing cosmic slrong wakes (Brandenberger el al. 
2010). 21cm observalions over wide sky areas offer a way of probing dislincl Hubble volumes in 
fhe same way as fhe CMB. Given fhe relalive infancy of Ihis field, if is likely lhal new ideas for 
probing cosmology will be developed during fhe developmenl of SKA. 

If is hopefully clear lhal SKA-LOW will have somelhing lo say on a wide range of cosmo¬ 
logical topics. The 21cm signal offers more varied routes to cosmology lhan Iradilional probes of 
large scale slruclure. This means lhal Ihe subllelies of exlracling precision cosmology from 21cm 
observations are nol fully understood and much work will be needed before SKA. SKA-LOW will 
Iransform our view of Ihe aslrophysics of Ihe epoch of reionizalion and cosmic dawn, Ihe same 
observations will add to our undersfanding of cosmology and provide new insighls into Ihe make 
up of Ihe early Universe. 
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